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摘 要 : 为 了 揭示 沙 区 建设 光伏 电站 后 沙 尘 运 移 特征 和 形态 发 育 及 演变 过 程 ,本 文 定 量 分 析 了 主要 风向 (光伏 板 
面向 正 南 ,风向 为 W, 即 光伏 设施 排 布 方向 与 风向 夹 角 为 0° ) 条 件 下 , 库 布 齐 沙漠 200 MWp 光伏 电站 腹地 区 域 光 伏 
板 不 同位 置 输 沙 特征 、 地 表 风 蚀 状况 和 流 场 分 布 规律 。 结 果 显示 :中 光伏 板 不 同 部 位 输 沙 率 随 风速 的 增 大 而 增加 ， 


不 同 风速 下 平均 输 沙 率 板 间 (1.17 g- cm? + min |) > 板 后 (0. 86 g cm? + min”! 
min!) O ~30 em 垂直 断面 上 ,不 同 部 位 90% 以 上 输 沙 率 集中 分 布 在 0 ~ 8 cm 高 度 范围 内 ,95% 以 上 输 沙 率 集中 


> 板 前 (0.65 g- cm + 


VS 


分 布 在 0 ~11 cm 高 度 范围 内 , 挟 沙 气流 中 输 沙 率 随 高 度 的 分 布 均 符合 指数 规律 递减 (RR 二 0.98)。@ 沙 粒 平均 路 


移 高 度 表 现 为 板 前 > 板 后 > 板 间 ， 


与 风速 正 相关 。 风 沙 流通 量 系 数 分 析 表 明 , 板 间 和 板 后 风沙 流 集中 在 近 地 层 ， 


而 板 前 风沙 流 有 向 高 层 移动 的 趋势 。@) 光伏 板 不 同位 置 摩 阻 风速 表现 为 板 前 (0. 562 0 ~0.596 0 m.s”!') > 板 后 
(0.331 2 ~0.436 0 m+ s~') > 板 间 (0.325 2 ~0.363 2 m .s”)。 而 且 板 前 干 沙 层 厚度 显著 高 于 其 他 位 置 ,土壤 可 
刨 性 增加 ,综合 作用 下 导致 板 前 发 生 强烈 掏 蚀 现象 ,观测 期 风蚀 深度 可 达 12.44 cm。 人 研究 可 为 沙漠 地 区 光伏 电站 


内 风沙 危害 的 科学 防治 提供 理论 支撑 。 
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风沙 流 是 指 含 有 沙 粒 的 运动 气流 , 挟 沙 气流 中 
所 搬运 的 沙子 在 搬运 层 内 随 高 度 的 分 布 特征 称 为 风 
沙 流 结构 避 。 研 究 风沙 流 结构 可 以 分 析 风 沙 地 瑶 
形态 的 形成 发 育 .演变 发 展 等 规律 ,在 认识 客观 规律 
基础 上 ,采取 有 效 的 措施 控制 或 促进 风沙 运动 ,才能 
达到 除 害 兴 利 或 转 害 为 利 的 目的 。 我 国 西北 地 
区 光照 强 ,辐射 充足 , 沙 区 土地 使 用 成 本 低 , 为 太阳 
能 光伏 产业 发 展 提供 了 条 件 ” 。 随 着 大 规模 太阳 
能 光伏 发 电站 的 兴建 ,施工 对 地 表 扰动 ,土壤 活化 为 
风沙 活动 提供 了 丰富 的 沙 源 ,而 光伏 设施 干扰 下 流 
场 格局 发 生变 异 ,打破 了 原 有 维持 地 表 平 衡 状态 的 
动力 分 布 ,地 表 风 沙 运动 规律 随 之 改变 。 有 研究 
表明 ,光伏 电站 边缘 区 气流 受到 电 板 倾斜 向 下 的 导 
流 作 用 , 集 流 加 速 并 形成 回旋 涡流 ,导致 光伏 板 下 出 
现 掏 蚀 现象 ,形成 以 光伏 板 下 沿 为 轴线 的 风蚀 坑 
( 沟 ), 当 气流 经 过 集 流 加 速 在 光伏 板 背 风 侧 扩散 后 
能 量 衰减 , 沙 粒 堆 积 在 板 间 形 成 堆积 沙 克 地貌 ,严重 
威胁 光伏 电站 安全 运营 。 目 前 ,在 光伏 产业 领 
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域 ,学 者 们 主要 对 光伏 发 电 组 件 材料 及 安装 方 
式 “” 光伏 组 件 电学 模型 与 仿真 ”” 、 降 沾 对 发 电 
效率 影响 ” \ 光 伏 面板 积 尘 清洁 '" 等 研究 较 多 ,而 
关于 沙 区 光伏 设施 干扰 下 地 表 侵 蚀 动力 过 程 . 风 蚀 
地 貌 形 成 演变 规律 等 方面 研究 相对 注 弱 。 
笔者 在 库 布 齐 沙漠 200 MWp 光伏 电站 试验 观 
测 中 发 现 ,光伏 电 站 内 风沙 运动 特征 的 观测 结果 显 
著 受 到 光伏 设施 排 布 方向 与 风向 夹 角 的 影响 。 根 据 
实地 观测 ,光伏 板 面向 正 南 东西 排 布 ,在 环境 风向 
为 偏 W 或 偏 E 风 条 件 下 ,密集 排 布 的 光伏 板 对 风向 
会 产生 “整流 作用 ” ,电站 腹地 板 间 近 地 表 则 表现 为 
W 或 上 风向 ,此 时 电站 内 光伏 板 周 围 侵 蚀 和 堆积 过 
程 强烈 ; 当 环 境 风向 为 偏 S 或 偏 N 风 时 ,光伏 阵列 
有 着 类 似 沙 障 作用 ,截流 阻 沙 效应 显著 ,电站 腹地 侵 
蚀 堆 积 过 程 表现 微弱 。 基 于 此 ,研究 光伏 设施 排 布 
方向 与 风向 不 同 夹 角 条 件 下 的 电站 内 近 地 表 流 场 变 
异 规律 和 输 沙 通 量 特征 ,有 助 于 进一步 认识 沙漠 地 
区 架设 光伏 板 后 地 表 形 态 变 化 的 发 育 过 程 和 演变 规 
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律 , 进 而 可 为 更 好 地 确定 光伏 基地 风沙 危害 防治 技 
术 方 案 提 供 理 论 依据 。 

本 研究 主要 定量 分 析 了 主要 风向 (光伏 板 面 向 
正 南 , 风 向 为 W, 即 光伏 设施 排 布 方向 与 风向 夹 角 
为 0") 条 件 下 , 库 布 齐 沙漠 200 MWp 光伏 电站 腹地 
光伏 板 不 同位 置 答 沙 通 量 AFL AVM EF J HOBO 小 
型 移动 气象 站 和 钢 尺 等 工具 观测 光伏 板 不 同位 置 的 
刨 积 状 况 .风速 变化 和 干 沙 层 厚度 等 指标 , 旨 在 揭示 
光伏 设施 排 布 方向 与 风向 夹 角 为 0" 条 件 下 电站 腹 
地 光伏 板 周 围 的 风沙 流 活动 规律 ,以 期 为 全 面 认 识 
沙 区 建设 光伏 电站 后 地 表 形 态 发 育 过 程 提供 依据 。 


1 研究 方法 


1.1 研究 区 概况 

研究 区 位 于 内 蒙古 鄂尔多斯 市 杭 锦 旗 独 贵 塔 拉 
镇 工业 园区 库 布 齐 沙漠 200 MWp 光伏 电站 内 ,地 理 
坐标 为 37"20' ~39°50'N ,107*10' ~111°45'E, ZK 
域 属于 温带 大 陆 性 气候 ,年 平均 气温 为 5~8 Y ,年 
降水 量 在 258.3 mm 左右 ,年 潜在 蒸发 量 2 400 mm 
左右 ,年 太阳 总 辐射 量 597.9 kJ . cm ,无 霜 期 149 
d, 风 沙 活动 主要 集中 在 3 一 5 月 ,全 年 大 风 日 数 为 
25 ~35 d, 全 年 盛行 西北 风 和 西风 ,夏季 盛行 东南 
风 。 平均 风速 二 5 m+ s 的 次 数 为 323.4 次 ,最 多 
年 份 达 418 次 。 该 地 区 的 沙丘 主要 呈 NW-ES 走向 
的 新 月 形 沙 丘 .新 月 形 沙丘 链 和 格 状 沙丘 链 等 ,沙丘 
高 度 10 ~60 m, 其 中 较为 平缓 的 沙 地 占 10% 左右 ， 
半 固 定 沙 地 占 30% ,60% 为 流动 性 沙 地 , 且 流 动 性 
较 强 。 

研究 区 光伏 板 设 置 方式 及 其 规格 :光伏 板 面向 
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下 南方 ,整体 东西 方向 排 布 , 相 邻 两 排 光伏 板 间 距 
900 cm, 面 板 上 沿 距 地 面 垂直 高 度 270 cm, 下 沿路 
地 面 垂直 高 度 35 cm。 单 组 光伏 板 由 2 排 18 列 99 
cm x195 cm 基本 光伏 板 单元 组 成 , 单 组 光伏 板 整体 
规格 为 400 cm x 1 800 cm, 地面 投影 宽度 约 为 320 
cm, 
1.2 试验 设计 与 数据 处 理 
1.2.1 试验 布设 开展 实验 前 光伏 板 不 同位 置 已 
经 发 生 严重 的 风蚀 和 堆积 现象 ,光伏 板 下 掏 蚀 形成 
以 光伏 电 板 前 沿 为 轴线 的 风蚀 沟 , 板 间 形成 堆积 沙 
EMEL) 。 因 此 ,在 观测 仪器 布设 前 ,首先 将 试 
验光 伏 板 周围 地 表 进 行 平整 , 整 平 后 的 光伏 板 下 沿 
高 度 距 离 地 表 约 35 cm( 图 2)。 

观测 仪器 具体 布设 方法 如 图 3 所 示 。 在 光伏 板 
板 间 、 板 前 和 板 后 3 个 位 置 同 时 放置 集 沙 仪 , 集 沙 仪 
采集 高 度 为 30 om, FE 15 个 进 沙 口 , 进 沙 口 规格 为 
2 cm x2 cm。 利 用 HOBO 小 型 移动 气象 站 同步 观 
测 光 伏 板 板 间 、 板 前 和 板 后 3 个 位 置 20 cm 50 cm, 
100 cm .200 cm 和 310 cm 高 度 处 风速 变化 情况 , 风 
速 数 据 记 录 间 隔 设 定 为 1 s, 数 据 采集 间隔 设置 为 
3 s。 同 时 利用 风蚀 测 钙 观 测 光 伏 板 周围 地 表 蚀 积 
状况 , 沿 垂直 光伏 板 方向 ,在 光伏 板 投影 面积 之 外 
( 即 光 伏 板 前 后 ) 测 钙 间 距离 100 cm ,进入 光伏 板 
投影 区 测 钙 间 上 距 为 SO cm; 这 样 的 布设 在 沿 平行 光 
伏 板 走向 布设 5 ÍT AL HW ET, EEFT M EF E E 200 
cm ; JET BE 80 cm ,地 上 外 露 40 cm ,地 下 埋设 40 


cm, 


试验 于 2019 年 3 月 25—27 日 进行 ,进行 了 风 
速 (风向 ) 、 输 沙 率 、 风 人 蚀 量 的 同步 观测 。 输 沙 观测 
依据 风 况 分 为 3 个 时 段 进行 ,每 次 观测 时 间 为 1 Ph， 


1 光伏 板 风蚀 地 表 
Fig.1 Wind erosion surface under the interference 


of PV panels 


图 2 ” 整 平 待 试验 光伏 板 地 表 


Fig.2 Smoothing underlying surface to test 
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Xi 面 图 900 cm 


风向 传感器 
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图 3 试验 布设 示意 图 
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Fig.3 Schematic diagram of the experimental layout 


观测 期 3 个 阶段 平均 风速 (2 m) X 7.79 m+ so! 
8.43 m.s 和 8.80m.s ,所 收集 到 的 沙 物质 用 
自封 袋 分 层 取样 , 带 回 实验 室 用 0.01 g 电子 天 平 称 
重 。 风 蚀 测 钙 测 定 为 一 场 风 ,时 间 段 为 3 月 25 H 
13:00 至 3 月 27 日 10:40。 
1.2.2 数据 处 理 

(1) 风沙 流 结构 特征 值 

为 阐释 单 组 光伏 板 不 同位 置 风沙 流 结构 特征 与 
沙 物质 吹 蚀 .搬运 和 堆积 的 关系 ,本 研究 引用 文献 
[11] 中 提出 的 风沙 流 结构 特征 值 入 ,作为 判断 地 表 
蚀 积 方向 的 指标 。 本 试验 采用 的 集 沙 仪 进 沙 口 为 
2 cm x2 cm, 大 量 研究 发 现 ,根据 拟 合 曲线 求 出 0 ~ 
1 cm 高 度 和 2 ~ 10 cm 高 度 的 输 沙 率 ,进而 推导 出 入 
值 并 判别 地 表 乌 积 方向 与 实际 观测 存在 的 偏差 。 本 
研究 根据 张 正 假 等 “研究 方法 ,利用 0 ~2 cm 和 
2 ~10 cm 高 度 输 沙 率 来 计算 其 风沙 流 特征 值 。 计 
算 公 式 如 下 : 


Q,.10 
r (1) 
式 中 :Oo 和 0 :分别 为 风沙 流 中 2 ~ 10 cm 高 度 
层 和 0 ~2 cm 高 度 层 的 输 沙 率 。 
(2) 赤 池 信息 量 准则 
风沙 流通 量 最 佳 拟 合 模型 采用 赤 池 信息 量 准则 


入 


(akaike information criterion ,简称 AIC) 确定 ,该 准则 
是 衡量 统计 模型 拟 合 效果 的 一 种 标准 ,AIC 值 越 小 ， 
表明 模型 拟 合 效果 越 好 。 计 算 公 式 为 : 


AIC =nln | È Cg; -4,)? | +2p (2) 


式 中 :n 为 集 沙 仪 集 沙 盒 个 数 ;9 为 集 沙 仪 实测 结 
果 ;4, 为 模型 拟 合 估算 结果 ;p 为 拟 合 模型 中 参数 的 
ABER 。 

(3) 净 风 蚀 量 计算 

根据 风蚀 测 镍 的 风蚀 与 沉积 计算 净 风 亿 量 417, 。 
XH, + EH, 


n pins 


IOP Ha 和 五 . 分别 表示 观测 点 的 沉积 和 风蚀 ,单位 
CM; 1 in BEAR LIN ST EFA PO" 。 

(4) 空气 动力 学 粗糙 度 和 摩 阻 速度 计算 

大 气 呈 中 性 或 近 中 性 稳定 条 件 下 ,地 表 风 速 廓 
线 一 般 满足 : 


Al, (3) 


U 
U,= yin (4) 


N 


式 中 ;U0, 为 摩 阻 速度 ;2 为 空气 动力 学 粗糙 度 ;K 为 
冯 卡 曼 常数 (0.4) ;2 为 距 地 面 高 度 0 。 
本 研究 采用 对 数 廊 线 拟 合法 计算 空气 动力 学 粗 


BEE , 测 得 5 个 高 度 (20 cm.、50 cm 100 cm 、 200 cm 

和 310 cm) 处 风速 ,用 最 小 二 乘 回 归 所 测 得 的 风速 

oy (5) 

式 中 :4.\B 为 回归 系数 。 在 公式 (4) 中 , 令 U, =0 可 
求 出 : 

Z, =exp( -A-B"') (6) 

由 公式 (3) 和 (5 ) 计算 可 得 摩 阻 速度 计算 方程 
Wi, 

U, =KB (7) 


2 结果 与 分 析 


2.1 光伏 板 不 同位 置 输 沙 率 特征 
通过 计算 光伏 板 不 同位 置 各 层 (0 ~30 em, f 2 
cm 一 层 , 共 15 层 ) 输 沙 率 ,然后 累加 得 出 总 输 沙 率 。 


结果 显示 ,不 同位 置 总 输 沙 率 随 风速 增 大 表现 出 增 
加 趋势 ,不 同 风 速 下 ,总 输 沙 率 均 表现 为 板 间 > 板 
后 > 板 前 。7. 79 m. s ”风速 下 板 间 (0. 953 g- 
cm .min-) > 板 后 (0. 570 ge cm + min ) > 
板 前 (0.496 g- cm .min-);8.43 m » s~! KUE F 
板 间 (1.074 g» cm” + min- ) > 板 后 (0.655 g 


cm” - min™') > 板 前 (0.602 g + cm + min- ); 
8.80 m + s7 风速 下 板 间 (1.470 g + cm~? + min™') 


> 板 后 (1.348 g + cm. min™') > 板 前 (0.847 g - 
cm + min”! )。 不 同 风速 下 平均 总 输 沙 率 为 板 间 
(1.17 g- cm .min-) > 板 后 (0. 86 g em + 
min`") > 板 前 (0.65 g .em + min ) , 板 前 和 板 后 
输 沙 率 分 别 为 板 间 的 55.62% 和 73.58% 。 

为 研究 输 沙 率 随 高 度 的 变化 规律 ,利用 Origin 
2017 做 出 散 点 图 (图 4) ,结合 前 人 的 研究 成 果 , 用 
指数 函数 9 = ae" AN FE PHBL q = aH" Fy DEFT TW 
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图 4 不 同 风速 下 光伏 板 不 同 部 位 输 沙 率 指数 函数 拟 合 


Fig. 4 Mass-ux-density proles at the different positions and velocity of PV panels 
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图 5 不 同 风速 下 光伏 板 不 同位 置 的 累积 输 沙 率 情况 


Fig.5 Cumulative percentage of sand transported rate at the different positions and velocity of PV panels 


合 , 式 中 9 为 输 沙 率 ,a 和 0， 为 系数 ,h 为 高 度 , 结 
发 现 不 同 风速 梯度 下 均 表现 为 指数 函数 相关 度 最 
高 ,R SO. 98。 说 明光 伏 板 不 同位 置 输 沙 率 随 高 度 
的 变化 更 加 符合 指数 函数 形式 , 拟 合 函 数 系数 用 来 
表征 不 同 条 件 下 沙 尘 通 量 的 变化 特点 。 对 于 指数 函 
数 而 言 , 有 研究 认为 系数 a 表征 风沙 流 中 沙 粒 浓度 
的 最 大 值 ,而 有 人 则 认为 是 蠕 移 输 沙 率 ;系数 5b 一 般 
认为 是 反映 沙 粒 浓度 随 高 度 的 衰减 程度 , 即 递减 率 ， 
b 值 越 大 衰减 速度 越 慢 ,b 值 越 小 衰减 速度 越 快 |。 
本 研究 系数 a 值 的 范围 为 0. 300 8 ~0.920 8 ,不 同 
风速 下 a 值 均 表 现 为 光伏 板 板 间 > 板 后 > 板 前 , 表 
明 板 间 近 地 表 风 沙 流 的 侍 粒 浓 度 或 蠕 移 量 大 于 板 后 
和 板 前 。 系 数 5b 值 范围 为 0. 279 3 ~ 0. 454 0 ,在 
7.97m…s 和 8.43 m-s 风速 下 系数 2 的 值 板 
前 > 板 间 > 板 后 ,表明 沙 粒 浓度 随 高 度 增加 递减 率 
依次 增 大 。 而 8.8 ms !' 风 速 下 系数 5 的 值 板 
间 > 板 后 > 板 前 ,此 时 板 前 沙 粒 浓度 随 高 度 变 化 递 
减 率 发 生 突变 。 
2.2 ”光伏 板 不 同位 置 沙 粒 的 跃 移 高 度 

为 更 加 直观 地 认识 光伏 板 不 同位 置 沙 粒 浓度 与 
高 度 之 间 的 关系 ,将 累积 输 沙 率 随 高 度 的 变化 做 成 
柱状 图 。 由 图 5 可 以 看 出 ,不 同 风速 条 件 下 沙 尘 输 
移 高 度 存在 差异 ,7.79 ms | ER PF 75% 以 上 的 
输 沙 率 集中 分 布 在 0 ~4 cm 高 度 , 除 板 间 位 置 以 外 ， 
90% 以 上 的 输 沙 率 分 布 在 0 ~6 cm 高 度 。 随 着 风速 
的 增加 ,不 同位 置 的 沙 粒 跃 移 高 度 随 之 升 高 ,8. 43 
ms-! 和 8.80 ms ! 风 速 下 75% 以 上 的 输 沙 率 在 
0 ~5 cm 高 度 范围 内 传输 ,90% 以 上 的 输 沙 率 分 布 
在 0 ~8 cm 高 度 。 比 较 不 同位 置 沙 粒 跃 移 高 度 差异 


性 ,7.79 m. s ”风速 下 总 体 表现 为 板 间 > 板 前 > HR 
后 ,8.43 m+ s 风速 下 板 前 沙 粒 跃 移 高 度 显著 上 


升 , 且 高 于 板 间 和 板 后 位 置 。 而 随 着 风速 的 继续 增 
大 , 板 前 沙 粒 的 跃 移 高 度 急 剧 下 降 , 低 于 板 间 和 板 后 


位 置 。 

平均 跃 移 高 度 是 表征 地 表 受 气流 剪 切 状况 的 参 
数 , 其 定义 为 在 该 高 度 以 上 及 以 下 气流 搬运 沙 物 质 
能 力 相同 , 即 累 积 输 沙 率 随 高 度 变 化 曲线 上 累积 输 
沙 率 为 50% 所 对 应 的 高 度 。 本 研究 中 7.97 mss! 
风速 条 件 下 沙 粒 平均 跃 移 高 度 板 前 (2.31 cm) > 板 
后 (2.25 cm) > 板 间 (1.98 cm) ;8.43 m © s~' RUBE 
条 件 下 沙 粒 平均 跃 移 高 度 板 前 (2.47 cm) > 板 后 
(2.22 cm) > 板 间 (2.08 cm); 而 8.80 m > s~! KUE 
条 件 下 沙 粒 平均 跃 移 高 度 板 间 (2. 26 cm) > 板 后 
(2.10 cm) > 板 前 (1.52 cm) 。 可 见 ,不 同 风速 条 件 
下 沙 粒 平均 跃 移 高 度 与 通 量 系数 反应 的 不 同位 置 沙 
粒 浓度 随 高 度 的 变化 趋势 相 一 致 ,在 7.97 m + 7 
和 8.43 m:s 风速 条 件 下 , 板 间 和 板 后 位 置 沙 粒 
平均 跃 移 高 度 低 ,系数 (EK b 值 小 ,风沙 流 运动 相 
对 更 加 贴近 地 表 。 而 板 前 沙 粒 平均 跃 移 高 度 高 , 系 
数 a 值 小 5 值 大 ,风沙 流 结构 有 向 上 层 移动 的 趋势 。 

综 上 所 述 ,在 8.80 m -s7 风速 下 , 板 前 风沙 流 
通 量 系数 5 EL 沙 粒 整体 和 平均 跃 移 高 度 均 出 现 突 
变 。 初 步 分 析 其 原因 ,是 随 着 观测 时 段 的 进行 ,在 
8.80 m+ s 风速 条 件 下 观测 时 , 板 前 地 表 已 经 发 生 
风蚀 现象 ,风蚀 深度 可 达 10 cm 以 上 , 而 地 表 风 蚀 导 
致 土壤 中 的 细 粒 物质 易 被 吹 蚀 ,相应 表层 粗 沙 组 分 
比例 增加 ,可 蚀 性 降低 。 从 输 沙 率 的 变化 中 也 可 
以 印证 这 一 事实 ,8.80 m + s ”风速 下 板 前 位 置 输 沙 
率 为 不 同 部 位 平均 输 沙 率 的 69.33% ,而 7.97 m - 
s Fil 8.43 m © s! 风速 下 分 别 为 73. 70% 和 
77.48% 。 同 时 地 表 粗 粒 化 ,使 得 沙 粒 平均 跃 移 高 度 
降低 ( 表 1) ,8.80 m + s7 风速 下 板 前 近 地 表 沙 粒 浓 
度 占 比 显著 增加 ,0 ~ 2 cm 层 输 沙 率 占 比 由 8. 43 
m:s 时 的 41. 73% 上 升 到 8.80 m > s ”时 的 
57.24% ,增加 了 15.51% ,而 近 地 表 沙 粒 浓度 占 比 
的 增加 导致 沙 粒 浓度 随 高 度 增加 递减 率 增 大 , 进而 


于 g 


表 1 光伏 板 不 同 风速 不 同位 置 0 ~2 cm 层 输 沙 率 
占 比 情况 对 比 
Tab.1 Proportion of sand transport rate at the 
height of 0 -2 cm at the different positions and 
velocity of PV panels 


风速 /( m» s~) 板 间 /% 板 前 /% 板 后 /% 
7.97 45.57 43.55 50. 12 
8.43 46.41 41.73 46.62 
8.80 45. 82 57.24 47.97 


出 现 系数 4b 值 急 剧变 大 , 沙 粒 整体 和 平均 跃 移 高 度 
均 变 低 的 现象 。 
2.3 ”光伏 板 不 同位 置 蚀 积 规律 

根据 风沙 流 特征 值 计算 公式 ,得 出 不 同 风速 下 
光伏 板 不 同位 置 风沙 流 结构 特征 值 。7.97 m+ s 
风速 条 件 下 , 板 间 、 板 前 和 板 后 的 入 值 分 别 为 1. 09、 
1.21 和 0.94;8.43 m+ s” 风速 条 件 下 , 板 间 、 板 前 
和 板 后 的 入 值 分 别 为 1.05.1.23 和 1.08; 在 8.80 
m- s ”风速 下 , 板 间 、 板 前 和 板 后 的 入 值 分 别 为 
1.05 .0.67 和 0.98。 综 合 来 看 ,光伏 板 板 间 的 入 值 
大 于 1 ,应 处 于 弱 风 蚀 状 态 。 板 后 的 入 值 接近 于 1, 
应 处 于 平衡 状态 。 板 前 除 8. 80 m- s ”风速 条 件 下 
和 值 大 于 1, 近 地 表 风沙 流 未 达到 饱和 流 ,造成 板 前 
地 表 强 烈风 蚀 。 板 前 8.80 m . s 风速 条 件 下 入 为 
0. 67 , 而 实际 观测 板 前 位 置 仍然 处 于 风蚀 搬运 状态 ， 
出 现 此 种 情况 的 原因 同样 由 于 风蚀 过 程 导致 地 表 粗 
BAL?) ,致使 蠕 移 量 加 大 ( 表 1) ,根据 A 计算 公式 
可 知 ,分 母 增 大 ,和 值 减 小 。 

利用 风蚀 测 钙 计 算 光 伏 板 不 同位 置 净 风 乌 量 发 
现 ( 图 6), 板 前 ( 即 光伏 板 下 沿 ) 附近 表现 为 强烈 风 
ih ,风蚀 深度 达 12.44 cm。 在 板 前 向 板 间 方向 距离 
光伏 板 下 沿 1.5 m 处 周围 表现 为 强烈 堆积 , 堆积 高 


光伏 板 下 沿 光伏 板 上 沿 
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图 6 垂直 光伏 板 方向 不 同位 置 净 风 人 蚀 量 差异 性 


Fig.6 Volume of erosion and deposition around 


the test solar PV panels 


度 达 7. 38 cm。 而 在 板 间 和 板 后 位 置 蚀 积 活动 较 
弱 ,这 与 风沙 流 结构 特征 值 分 析 结果 基本 一 致 。 
因此 ,光伏 设施 排列 方向 与 风向 夹 角 为 0" 情况 下 板 
前 会 出 现 强烈 掏 蚀 现象 ,而 板 间 和 板 后 的 地 表 蚀 积 
变化 微弱 ,可 见 , 板 前 位 置 是 沙 区 光伏 电站 内 防治 的 
重要 部 位 。 


3 讨论 
3.1 风沙 流通 量 


不 同 下 执 面 近 地 表 风沙 流通 量 究 竞 符合 怎样 的 
函数 分 布 至 今 仍然 没有 定论 。 刘 芳 等 "对 乌 兰 布 
和 沙漠 东北 缘 5 种 典型 下 热 面 0 ~ 100 cm 内 输 沙 率 
随 高 度 的 变化 进行 函数 拟 合 , 发 现 短 函 数 的 相关 性 
最 好 。 吴 晓 旭 等 "对 毛乌素 沙 地 风沙 流 结构 研究 
发 现 ,平坦 流沙 地 ,流动 沙丘 与 固定 沙丘 丘 顶 上 输 沙 
率 随 高 度 都 呈 短 函数 分 布 , 半 固 定 沙丘 丘 顶 的 输 沙 
率 随 高 度 呈 指 数 分 布 。 张 正 候 等 "对 腾 格 里 沙漠 
东南 缘 风 沙 流通 量 进 行 了 长 期 观测 ,利用 5 种 常用 
拟 合 模型 对 比 研究 发 现 , 平 坦 沙 地 上 风沙 流通 量 更 
符合 指数 函数 分 布 。 陈 曦 等 * 对 乌 兰 布 和 沙漠 东 
南 缘 光 伏 电 场 电 板 的 前 沿 . 后 沿 及 光伏 阵 行道 处 风 
沙 流 结构 研究 表明 , 电 板 不 同 部 位 输 沙 率 随 高 度 增 
加 符合 多 项 式 函 数 形式 。 本 研究 基于 野外 观测 数 
据 , 采 用 指数 函数 和 短 函 数 进行 拟 合 分 析 , 运 用 赤 池 
信息 量 准则 研究 发 现 ,指数 函数 的 相关 度 最 高 , 且 
AIC 值 显著 低 于 其 他 函数 模型 ( 表 2)。 因 此 ,指数 
函数 能 够 很 好 模拟 库 布 齐 沙漠 200 MWp 光伏 电站 
内 风沙 流通 量 。 

观测 方法 的 不 同 所 得 出 系数 a.b 值 的 范围 不 
同 。 对 于 指数 函数 而 言 , 系 数 5 值 受 拟 合 模型 高 度 
单位 的 影响 ,利用 高 度 单位 cm 拟 合 时 , 相 较 于 利用 
高 度 单位 m 进行 拟 合 , 系 数 b 值 扩大 100 倍 。 函 数 
式 选择 的 不 同 也 是 影响 4 值 大 小 的 因素 ,目前 用 的 
站 数 函 数 模 型 有 q(h) =abi q(h) =ae q(h) = 
。 计 算 可 得 5 =e" =e" ,可 见 , 讨 论 系 
数 b 值 范围 时 应 在 统一 函数 式 的 基础 上 。Dong 
等 .通过 风 洞 试 验 研 究 表明 ,系数 a 的 值 小 于 2. 5, 
b 值 最 大 可 达 0. 35 [高 度 单位 为 em,g (h) = 
ae“* "]。 哈 斯 中 通 过 野外 实验 观测 得 出 ,4 ~ 10 
m 高 的 横向 沙丘 系数 a 值 小 于 13,[0,1] 区 间 值 占 
总 体 的 70% ,b 值 范 围 为 0.193 ~0.494[ 高 度 单位 
为 cm,g(h) =a”), REJ b 值 最 大 值 为 5. 18 
[g(h) =ae** ']。 张 正 候 等 在 平坦 沙 地 通过 


ae =h. by! 
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表 2 输 沙 率 随 高 度 变化 的 指数 函数 和 窘 浮 数 回归 结果 对 比 
Tab.2 Results of the exponential and power regression analyses for the mass-ux-density profiles 
风速 观测 指数 函数 拟 合 FE PRU 
/( m+s) 位 置 公式 R AIC 公式 R AIC 
7.97 板 间 Q =0.596 Se ~% 39 8 0.998 5 -117. 11 Q =0. 212 5h~°™? 0.5179 -5.88 
板 前 Q =0.300 8e "2 8 0.983 8 - 100. 42 Q =0. 110 9h~°%? 0.516 5 一 24. 95 
板 后 Q =0.411 le-0348 3 0. 995 9 -114. 47 Q =0. 130 1h-°%! 0.465 5 -20. 68 
8. 43 板 间 Q =0. 682 8e 70312 4h 0. 999 9 - 153. 36 Q =0. 520 6h 7% 867 0. 935 9 -56. 68 
板 前 Q =0. 342 2e 70-77 3A 0.991 2 - 105. 60 Q =0. 268 1h~°°'88 0. 893 8 - 68. 28 
板 后 Q =0. 428 4e 70-313 8 0. 993 0 -103.719 Q =0. 322 3h °°? 0. 906 2 -64. 72 
8. 80 板 间 Q =0. 920 8e 7°31 2 0. 999 7 - 129. 49 Q =0. 702 1h~° °° 0.941 8 -49. 35 
板 前 Q =0.757 3e 70454 0h 0. 997 3 — 106. 75 Q =0. 492 4h71 184! 0.980 5 -77.16 
板 后 Q =0. 908 2e ~% 3322 0.999 1 - 113. 44 Q =0. 674 0h 7+9 0. 937 2 一 49. 28 


长 期 野外 实验 观测 发 现 ,系数 v 值 分 布 范围 较 大 ,最 
大 值 可 达 547 左右 ,系数 b 值 最 大 值 为 0.06[ 高 度 
单位 为 m,g(h) =ae““”] , 即 em 高 度 单位 下 拟 合 
最 大 值 可 达 6。 本 研究 试验 下 垫 面 为 无 植被 覆盖 的 
流动 沙 地 ,光伏 设施 干扰 下 不 同位 置 风沙 流通 量 系 
数 a 值 小 于 0. 92 ,系数 值 在 2.2 ~3.6 之 间 [ 高 度 
单位 为 cm,q(h) =ae“““ ] ,系数 a 值 的 范围 与 哈 
H Dong 等 .研究 结果 基本 一 致 ,系数 b 值 的 范 
围 略 小 于 哈 斯 55 、 张 正 假 等 "的 野外 观测 结果 。 

风沙 运动 是 一 种 贴近 地 表 的 沙子 搬运 现象 。 本 
文 基于 野外 观测 发 现 ,不 同 风速 下 光伏 板 不 同位 置 
75% 以 上 的 输 沙 率 集中 分 布 在 0 ~5 cm 高 度 ,90% 
以 上 的 输 沙 率 集中 分 布 在 0 ~8 cm 高 度 ,95% 以 上 
的 输 沙 率 集 中 分 布 在 0 ~ 11 em 高 度 。Chepil'” 发 
现 ;土壤 表面 0 ~ 5 cm 高 度 内 输 沙 率 占 总 体 的 
60% ~80% ,这 与 本 研究 结果 相 一 致 。 吴 正 “ 野外 
观测 表明 ,气流 挟 带 沙 量 90% 以 上 在 地 表 30 cm 高 
度 内 输送 ,其 中 约 有 80% 输 沙 率 在 0 ~ 10 cm 内 传 
o REKO 通过 实验 研究 认为 ,90% 的 输 沙 率 在 
0 ~10 cm 高 度 内 传输 。 张 正 候 等 ”通过 451 次 风 
沙 流 观 测 数据 研究 得 出 , 腾 格 里 沙漠 东南 缘 平坦 沙 
地 上 0 ~20 cm 高 度 范围 内 输 沙 率 占 总 输 沙 率 的 比 
例 大 于 95% 。 陈 新 闻 等 对 乌 兰 布 和 沙漠 流动 沙 
丘 时 外 观测 得 出 ,80% 以 上 输 沙 率 集 中 于 0 ~ 10 em 
高 度 ,90% 以 上 的 输 沙 率 集中 于 0 ~20 cm 高 度 。 这 
些 研究 结果 均 略 微 高 于 本 研究 ,可见 ,光伏 设施 的 存 
在 ,使 得 电站 内 风沙 流 垂 直 结 构 发 生变 化 ,与 裸 沙 地 
相 比 有 向 下 层 移动 的 趋势 。 
3.2 ”光伏 板 不 同位 置 蚀 积 机 理 

通过 计算 光伏 板 不 同位 置 风沙 流 结构 特征 值 发 


现 , 板 间 处 于 弱 风 蚀 状 态 , 板 后 处 于 弱 堆 积 状 态 , 板 
前 处 于 强风 蚀 状 态 ,而 且 风 蚀 测 钙 观 测 数据 也 很 好 
的 佐证 了 风沙 流 结构 特征 值 所 得 出 的 结果 。 引 起 光 
伏 板 不 同位 置地 表 蚀 积 状态 改变 的 主要 原因 是 ,由 
于 光伏 板 的 存在 改变 了 近 地 表 流 场 格局 ,已 有 研究 
表明 边缘 区 光伏 电 板 干扰 下 形成 板 下 集 流 加 速 区 ， 
板 前 板 后 阻挡 减速 区 , 板 面 抬升 加 速 区 , 板 间 消散 恢 
复 区 ,而 腹地 区 域 光 伏 电 板 强烈 影响 两 侧 流畅 格 
局 “” 。 本 研究 对 库 布 齐 沙漠 200 MWp 光伏 电站 
腹地 区 域 风 速 流 场 进行 了 观测 ,结果 显示 不 同 风速 
条 件 下 光伏 板 板 间 、 板 前 、 板 后 风速 随 高 度 增加 变化 
规律 趋 于 一 致 , 呈 " 卫 形变 化 (图 7) ,其 中 板 前 受 》 
伏 组 件 影响 ,风速 随 高 度 变化 幅度 较 大 , 板 间 和 板 后 
20 ~ 100 cm 高 度 范围 内 风速 变化 急剧 ,而 100 ~310 
cm 高 度 范围 内 风速 变化 平缓 。 可 见 ,光伏 板 的 存在 
主要 影响 近 地 表 100 cm 高 度 范 围 内 的 风速 流 场 ， 
100 cm 高 度 范围 内 风速 大 小 均 表 现 为 板 间 > 板 
后 > 板 前 ,这 与 不 同位 置 总 输 沙 率 大 小 关系 相 一 致 。 
计算 光伏 板 不 同位 置 空气 动力 学 粗糙 度 发 现 , 板 前 
(0.848 2 ~1.1000 cm) > 板 后 (0.036 1 ~0. 086 9 
cm) > 板 间 (0. 008 8 ~0.018 2 cm) ,这 与 近 地 表 
100 cm 高 度 范 围 内 不 同位 置 风速 大 小 呈 反 比 关系 ， 
与 不 同位 置 总 输 沙 率 呈 反比 关系 ,说 明光 伏 板 主要 
影响 了 板 前 周围 气流 场 , 增 加 地 表 粗 糙 度 的 同时 减 
小 了 输 沙 率 。 

计算 光伏 板 不 同位 置 摩 阻 风速 表现 为 板 前 
(0.562 0 ~ 0. 596 O0 m » s~!) > 板 后 (0.331 2 ~ 
0.436 0 m-s'') > 板 间 (0.325 2 ~ 0.363 2 m- 
s”) 。 表 明 板 前 地 表 所 受气 流 剪 切 力 最 大 ,这 与 平 
均 路 移 高 度 所 表征 的 结果 相 一致 。 同 时 , 本 研究 观 
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图 7 不 同 风 速 下 不 同位 置 风速 廓 线 


Fig.7 Wind velocity proles at the different positions and velocity of PV panels 


测 了 试验 期 间 光 伏 板 不 同位 置 的 干 沙 层 厚度 ,发 现 
板 间 位 置 平 均 为 10. 17 cm ,标准 差 0. 65 。 板 前 位 置 
平均 为 25.57 om, 标准 差 1. 00。 板 后 位 置 平均 为 
13.73 ecm, 标准 差 2.07。 干 沙 层 厚度 在 不 同位 置 差 
异性 较 大 , 板 前 和 板 间 干 沙 层 厚度 差异 达到 15 cm 
以 上 。 干 砂 层 厚 度 的 不 同 表 明光 伏 板 不 同位 置 土壤 
表层 含水 量 存在 较 大 差异 ,土壤 含水 量 达 到 一 定 阔 
(EJE , 沙 粒 间 的 黏 结 力 增 大 ,起 动 风 速 也 随 之 增 大 ， 
Sach UNO! 。 综 上 所 述 ,光伏 板 干扰 下 板 前 
地 表 受 气流 剪 切 力 增 大 和 表层 土壤 可 亿 性 增加 可 能 
是 引起 地 表 强 烈 掏 蚀 的 主要 原因 。 


4 结论 


(1) 库 布 齐 沙漠 200 MWp 光伏 电站 腹地 区 域 
光伏 板 不 同位 秆 输 沙 率 大 小 关系 为 板 间 > 板 后 > 板 
前 ,0 ~30 cm 高 度 范 围 内 挟 沙 气流 中 输 沙 率 随 高 度 
的 分 布 均 符合 指数 规律 递减 (RR SO. 98 ) 。 通 过 对 
光伏 电站 腹地 区 域 光 伏 板 不 同位 置 沙 粒 跃 移 高 度 和 
通 量 系 数 "综合 分 析 得 出 , 板 间 风沙 流 运动 相对 
更 加 贴近 地 表 ,而 板 前 有 向 上 层 移动 的 趋势 。 

(2) 结合 不 同位 置 风 况 观 测 结果 , 板 前 位 置 受 
到 光伏 设施 的 阻 滞 作 用 最 为 强烈 , 近 地 表 100 cm 高 
度 范 围 风速 明显 低 于 板 后 和 板 间 位 置 ,但 是 摩 阻 速 
度 增 大 , 且 沙 粒 的 整体 和 平均 跃 移 高 度 显 著 高 于 板 
间 和 板 后 ,表明 板 前 位 置 气流 对 地 表 前 切 力 较 大 , 同 
时 光伏 板 干扰 下 板 前 位 置 干 沙 层 厚度 增加 ,土壤 表 
层 含水 量 下 降 ,可 蚀 性 上 升 , 从 而 导致 了 板 前 位 置 处 
于 强风 蚀 状态 。 风 沙 流 结构 特征 值 和 风蚀 测 钙 的 实 
际 观测 数据 结果 也 很 好 地 佐证 了 这 一 结果 。 

(3) 本 研究 初步 证 实 光 伏 电 站 腹地 区 域 和 边缘 
区 域 侵蚀 机 理 异 同 ,但 关于 光伏 板 干 扰 下 的 侵蚀 动 


力学 机 制 任 然 不 能 完全 解释 ,有 必要 继续 开展 光伏 
设施 干扰 下 流 场 格局 变异 及 其 与 地 表 侵 蚀 堆积 过 程 
动力 关联 的 系统 研究 ,明晰 地 表 形 态 演化 的 过 程 ,可 
为 沙漠 地 区 光伏 电站 内 风沙 危害 的 科学 防治 提供 理 
论文 撑 。 
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facilities in sandy areas 


TANG Guo-dong', MENG Zhong-ju', GAO Yong', ZHENG Hai-feng’, JIA Rui-ting', SHI Tao! 
(1. Desert Science and Engineering college/Key Laboratory of Aeolian Sand Physics and Sand Control Engineering 
in Inner Mongolia, Inner Mongolia Agricultural University , Hohhot 010011 , Inner Mongolia , China; 
2. Shenhua Group Zhungeer Energy Co. ,Ltd , Ordos 010300 , Inner Mongolia , China) 


Abstract; This study assesses the characteristics of sand transport, morphological development and their evolution 
after the construction of a photovoltaic (PV) power station in a sandy area. Field wind flow,sediment transport and 
surface elevation changes were analyzed in the hinterland of a 200 MWp PV power station in the Hobq Desert, 
where the angle between the arrangement of the PV facilities and wind direction is 0° (i. e. ,the PV panels face S 
and the wind direction is W). The results show that; (1) The total sand transport rate sequences around the PV 
panels were highest in the zones between panels (1.17 g + cm7? + min!) , followed by the zones behind panels 


(0.86 g +> cm 


=) 


- min~') ,and finally the zones in front of panels (0.65 g + cm~? + min™'). Sand transport rate 
was also found to decrease with increasing height. On the vertical section of 0 -30 cm ,more than 90% of the total 
sand transport flux was distributed within the height 0 -8 cm,and over 95% of the total sand transport flux was dis- 
tributed within the height 0 ~ 11 cm. The sediment discharge over three typical photovoltaic panel positions de- 
creased with height and conformed to an exponential function distribution (R? 0.98). (2) The average saltation 
height was positively correlated with wind speed; saltation height was, in order from highest to lowest, in front of 
panels > behind panels > between panels. The sand flux coefficient of wind-sand flow was also analyzed, and it was 
found that sand flux tends to be highest in front of panels , whereas between and behind panels, the flux was mostly 
concentrated in the near ground layer. (3) Under the interference of PV panels, shear velocities were 0.562 0 - 
0.596 0,0.331 2 -0.436 0 m+ s~',and 0. 325 2 -0.363 2 m+ s™' for the zones in front of ,behind ,and between 
panels , respectively. The thickness of the dry sand layer of the zones in front of panels was also significantly higher 
than at other locations , increasing soil erodibility. A strong erosion occurred in the zones in front of panels by com- 
prehensive action ,and wind erosion depth reached 12. 44 cm during the observation period. This study provides a 
better technical scheme for wind-sand hazards for the case of solar PV power stations in order to ensure their stable 
and safe operation. 

Key words; sand transport; aeolian sediment flux structure; fitted model; wind profile; photovoltaic power sta- 


tion; the Hobq Desert 


